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In the introduction, we will discuss the symmetries of QCD, its breaking in the physical 
vacuum and its restoration at finite temperature and/or density. Then we will discuss 
what exact relations can be made between chiral symmetry restoration and the imaginary 
part of the vector vector correlation function, which is directly related to the dilepton 
spectrum coming from a hot and/or dense medium. Finally, we will focus on the vector 
meson pole appearing in the imaginary part and discuss results on how the shape changes 
using QCD sum rules in medium. In particular, we will discuss three characteristic changes 
of the vector meson peaks; namely, shift of the peak position, increase of its width and 
Qh! the dispersion effect due to the breaking of Lorentz invariance in medium. 



1. INTRODUCTION 

At zero temperature and/or baryon density, the chiral symmetry is spontaneously bro- 
ken (qq) ^ and the hadronic world has only SU(N)y x U(l)y symmetry. However, if 
the temperature or baryon density increases, there will be a critical boundary above which 
the broken chiral symmetry gets restored (qq) = and the hadronic matter changes to a 
quark gluon plasma (QGP) state. This QGP state will have SU(N) L x SU(N) R x U(l) v 
symmetry, which is also the symmetry of QCD at the quantum level. 

1.1. Symmetry of QCD 

The Lagrangian of QCD is; 

Cqcd = ~-F; v F; v - £ q f (ip + M f )q f (1) 

where F* v = d^Al — d u A a ^ + igf abc 'A h ^A c v and A a ^ represents the spin 1 gauge field with 
color index a. qf represents the spin | matter field which has left and right chirality 
components, qi,R = §(1 =F 75)?, and the index / represents the N number of flavors such 
as the u,d,s- -. In the chiral limit (M — > 0), Cqcd is symmetric under U(N) L x U(N) R , 
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i.e., the lagrangian is symmetric under the following left right transformation. 
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where Q = —^-(Lq + LiX 1 ) + -4p-(i2 + -Rj-V). Under this transformation, the chirality 
does not mix (q L)R -> q L)R ). 

However, at the quantum level, due to the axial anomaly, the Ua(1) symmetry (Q = 
75 {Ro ~ Lq)) is broken and QCD is symmetric only under 

SU(N) L x SU(N) R x U(l) v . (3) 
This is also the symmetry of the QGP state. 

1.2. Spontaneously broken chiral symmetry 

In the hadronic world, the QCD symmetry at quantum level is further broken down 
spontaneously to the following group. 

SU(N) V x U(l) v (4) 

As a consequence the QCD vacuum and physical excitations are not symmetric under 
Q = l§Q a ^ a transformation. This is manifest from looking at various order parameters 
of chiral symmetry. 

First, consider the vacuum state. The chiral partner of any quark composite oper- 
ator Op is obtained by taking its commutator with the chiral charge operator F§ = 
J d 3 xq r yo'~f 5 ^-q. Hence the chiral partner of qq is q'Jz^-q. Since the vacuum expectation 
value of the latter vanishes identically, (qq) ^ shows that chiral symmetry is broken in 
the vacuum. The fact that chiral symmetry is spontaneously broken is also manifest in 
the presence of massless Nambu-Goldstone Bosons ( the 7r's), which couple to the broken 
axial current (0|A^|7r) ~ /„.. 

As can be seen from table 1, the broken chiral symmetry is also manifest in the mesonic 
excitation, such as m p ^ m ai , and baryonic excitations, such as mp ^ ^SulUJH- 



Table 1 



Chiral order parameters and evidence for its breaking in hadronic world. Here V£ = 
qjn^-q, A® = qj^^-q. ^_ is the chiral partner of \& + . 



Order parameter 




Hadronic world 


(qq) 




(QQ) ^ 0,/^ m 2 = _( m gg) 


W W(o)> - 


(A«(x)A«(0)) 


m p (770) ^m ai (1250) 


(* + (x)tt+(0)) " 


(*_(ar)*_(0)) 


mp (938) ^ m 5ll (1535) 
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1.3. Chiral symmetry restoration at high temperature and/or density 

As seen from lattice gauge theory calculations, the broken chiral symmetry will be 
restored at high temperature. The transition will become smooth as we increase the 
current quark mass. The actual order of the phase transition at the physical value of the 
quark masses for the u, d, s quarks are still controversial. The result using the Staggered 
fermions|3| shows a smooth transition at the physical value of the quark masses. However, 
the result using the Wilson fermions[Q shows a first order transition. 

QGP 

Chiral symmetry is restored 




5-10 n n /n 



Figure 1. Phase diagram of QCD. 

At present, the lattice calculation at finite density suffers a technical difficulty due to 
complex fermion determinant. Nevertheless, many model calculations, such as instanton 
models|5|], NJL model, Random Matrix model|| suggest that chiral symmetry will also 
be restored at finite baryon density. In fact, already at nuclear matter density, there exist 
model independent relations that shows that the chiral order parameter will be reduced to 
about 80 % of its vacuum value. Therefore, any physical effect related to chiral symmetry 
restoration will also be partially visible in nuclear matter. 

1.4. ^a(I) restoration at high temperature and/or density 

Present lattice calculations shows that the £7a(1) symmetry is not restored above the 
chiral phase transition J7|||]. This is reasonable because the {/a(1) breaking effect comes 
from a totally different mechanism from that of chiral symmetry breaking. Nevertheless, 
the present lattice calculations were performed with N = 2 flavors and studied the £/a(1) 
order parameter constructed out of two quark bilinears (2-point function). It would be 
interesting to perform similar calculations with larger number of flavors. The reason, 
which is valid above T c , is the following. The Ua{\) breaking effect comes from topo- 
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logically non-trivial gauge field configurations. Its contribution to the partition function 
has measure zero in the chiral limit, i.e. its contribution is proportional to the M N p inst , 
where N is the number of quark flavors having current quark mass M, and pi nst is the 
pseudo-instanton density||. For U^(l) order parameter, constructed out of n quark bi- 
linears, the breaking contribution is proportional to M N ~ n p inst . This implies, that in the 
chiral limit, after chiral symmetry restoration, the C/a(1) breaking effect shows up only 
when n > AT pp. 



Table 2 



?7a(1) breaking operators and its breaking contribution from the partition function. 
Here n a = qjs^q, S a = q^q. N is the number of quark flavors with mass M. p inst is the 
pseudo-instanton density. 



Breaking operator 


Breaking contribution 


(qq) 


M N ~ l Pinst 


(n a (x)n a (0)) - (S a (x)5 a (0)) 


M N ~ 2 p inst 


n-point function 


M N ~ n p inst 



2. OBSERVING CHIRAL SYMMETRY RESTORATION 

Among the order parameters of chiral symmetry restoration, such as those discussed 
in Table. 1, let us concentrate on the difference between the vector and axial vector two 
point functions. Using the dispersion relation, it can also be related to the difference 
between the imaginary parts. 

where, 

Il£ = J d'xe^ (^Or)K(O)) , U% = J d A xe^(A,{x)A^)). (6) 

The imaginary part of the vector correlator can be observed from the hadronic part of 
the cross-section e + + e~ — > even 7r Fig.2 (a). This has a sharp rho meson peak (770 
MeV) followed by resonances, which add up like a continuum. The imaginary part of the 
axial correlator can be observed from the hadronic part of the decay r — > odd 7r Fig.2 
(b). This part has the broad ai peak (1260 MeV). 

The two imaginary parts are clearly different and is one of the experimental signature 
that chiral symmetry is broken in the vacuum. 

In the vacuum, due to current conservation and Lorentz invariance, the polarization 
functions in Eq. @ has only one independent function, which is a function of the four 
momentum q 2 = u 2 — q 2 . At finite temperature or density, the Lorentz invariance is 
broken and the polarization functions will have a transverse (IIt) and a longitudinal 
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Figure 2. a)e + +e — > even ir, b)r — > odd 7r. 



(IIl) component (with respect to q), which will depend on both the energy u and three 
momentum q. 



n M ^,q) = n T g 2 p^ + n ig 2 p 



(7) 
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where for q = (uj, q) and medium at rest, we have, 

- qiOj/q 2 , and P^ = (q^qy/q 2 - flf^ - PjJ . n L = n T only when q = 0. 

Above the critical temperature or density, chiral symmetry is restored and Imn^ T = 
Imllf T , which at sufficiently high temperature or density can be calculated using pertur- 
bative QCDjllJ|. At high temperature or density, the poles disappear and become part 
of the continuum. Now the question is how would the poles behave and approach the 
continuum from below the critical boundary? 

In fact, one can observe Imll y from a hot or dense matter by looking at the inverse 
process e + + e~ — > even ir. That is, by looking at the dileptons coming out from the 
medium^f. Since, the dileptons interact weakly in the strong interaction time scale, 
once they are produced, they will carry the information about Imll^ directly to detectors. 
This would be true in the time scales involved in Relativistic Heavy Ion collisions or 
experiments of vector meson production inside a Heavy Nucleus planned at GSI and 
KEK. Unfortunately, the inverse process of r — > odd 7T, that is r + v coming out from the 
medium is not experimentally observable and one can not directly observe ImH" 4 . So the 
best we can do to observe chiral symmetry restoration (at least in heavy ion collisions) is 
to experimentally observe Imn y (especially the characteristic changes in the peak) and 
do QCD model calculations to find its relation to chiral symmetry restoration. 

The dilepton production rate from a hadronic system at finite temperature is given as 
follows. 



dR 



2ai 



dM 2 3vr 2 



dk 



uj M 2 



[Imn L (cu, k) + 2Imn T (cu, k)}n{uj) 



(8) 
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where, u 2 = k 2 + M 2 and n(ui) is the thermal Boltzmann factor. Here, IraQ^r^, k) is the 
hadronic part of the electro-magnetic vector current correlation function. If the system 
evolves in a RHIC, it should be folded with the evolution of the system. 

Assuming no change in the hadron properties at finite temperature or density, the 
expected dilepton spectrum from a hot dense system will have a peak near the vector 
mesons. Now, if the vector mesons change their properties in medium, changes will occur 
in the peak. The three characteristic changes we will discuss here are, i) peak (mass) 
shift, ii) increased decay width, iii) dispersion effect (or finite q) effect. 

2.1. Model calculation of ImII y at finite temperature and density 

Many model calculations have been performed to calculate the changes of the vector 
meson properties at finite temperature or density. The models are either based on ef- 
fective quark lagrangian or the effective hadron lagrangian. The former approach have 
model dependence related to modeling the non-perturbative nature of QCD into the quark 
language. The latter have large model dependence in extrapolating any effective hadron 
lagrangian to the on shell point of the vector meson mass region. Table 3 summarizes the 
results of the model calculations. 



Table 3 

Model calculation for Imlly at finite temperature and density. The downward, upward 
and horizontal arrows show respectively the decreasing , increasing and little change of 

values. 

model mass width dispersion 



Walecka model fT3" | 



Quark- meson coupling model ||14 



Finite 
density 



Phenomenological meson exchange |T5 



Brown- Rho scaling ||16 



7T 



N — A dynamics JT7 



V — N — N* dynamics |8 



Chiral model with Vector mesons [19,|20] 



Chiral model with Low T theorem||21 



Finite 

temperature 



q-model with VMD]22l 



Hidden local symmetry p3 



Collision broadening p 



Massive YM with baryons||25 



2.2. QCD constraints on Imn 

Let us now try to construct QCD constraints on Imn v emphasizing its connection to 
chiral symmetry restoration. Consider either the transverse or the longitudinal polariza- 
tion function in Eq.(|7]). The starting point is the energy dispersion relation at finite q. 
For small q 2 < u 2 , we can make a Taylor expansion of the correlation function, which is 
even function of u and q, such that, 

ReIl L , T {u, q) = Re (ll {u, 0) + Ul tT (u, 0) q 2 + ••) 
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(9) 



where we have assumed the retarded correlation function. 

The real part of eq.(|) is calculated via the Operator Product Expansion (OPE) at 
large — uj 2 — > 00 with finite q. The full polarization tensor will have the following form. 



n 



[IV 



,^q) 



n.m. 



+ E 7^-2 hrg^ 1 ■ -q^A^in.m.) + ••] , (10) 

s,t=2 y 

where, Q 2 = q 2 — uj 2 . Here, A d,T {n.m.) represents the nuclear matter expectation value 
of an operator of dimension d and twist r = d — s, where s is the number of spin index. 
These operators are defined at the scale Q 2 and the c's are the dimensionless Wilson 
coefficients with the running coupling constant. This way of including the density effect 
is consistent at low energy f26|j. The first set of terms in eq.(|T0"D come from the OPE of 
scalar operators, the second set from operators with non zero spin s. From this general 
expression, we can extract the OPE of II , n]^ T . Including the contribution up to dim.6 
operators only, they will have the following form, 

n°H = n pert (u) + —X + —b° 6 

where ^(frg) in II are contributions from dim.4 (dim.6) operators and feg(&g) in II 1 are 
from dim.4 (dim.6) operators. It should be noted that in calculating II 1 , it is enough to 
calculate only the "non-trivial" q 2 dependence ; i.e. neglect the dependence coming from 
expanding Q 2 = q 2 — uj 2 . That is so because we are interested only in the Lorentz breaking 
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effect and not in the trivial q 2 dependence ||27|]. Eq.(|TT|) are asymptotic expansions in 1/ 
Looking back at the dispersion relation in Eq.(9), we can look at different moments of 
the imaginary part and equate it to each term in the OPE in Eq.(Pl). 



j du 2 (Imn(u) - \mR pert (u)) u n ~ 2 = -b n (12) 

Here, Imn(tt) per( is the corresponding imaginary part of H per t- For Imn°, we have 3 
constraints from n = 2,4,6 with 62 = 0, and for Imn 1 , we have 3 constraints from n = 
4, 6, 8 with 64 = 0, Imn pert = for both the longitudinal and transverse directions. These 
constraints are model independent constraints and relates QCD perturbative calculations 
and expectation values of local operators to the phenomenologically observable spectral 
density Imn(w). The usefulness of the constraints can be summarized as follows. 

1. The temperature dependence (density dependence) of the local condensates (b's) 
can be calculated using lattice QCD. Therefore, starting from the vacuum form of 
Imn(u), one can study its changes at finite temperature or density. 
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2. One can look at the effect of each local operators in b's separately. Hence, from 
looking at the effect of chiral symmetry breaking operators, one can study the effect 
of chiral symmetry restoration on Imll('u). 

3. All model calculations of ImII(w), or even the experimental data itself should satisfy 
the above constraints. 

It is straightforward to extract the chiral symmetry breaking operators. First, it should 
be noted that all the twist- 2 and Gluonic operators, appearing in Eq.(ll), are chirally 
symmetric operators. Four quark operators have a chirally symmetric component and a 
breaking part. Assuming a four quark operator O, one can divide these parts by subtract- 
ing and adding its chiral partner O = \[0 — (9(chiral partner)] + + O (chiral partner)]. 
The explicit form of b's can be found in [^j for Imll and in |27j for Imll 1 . It would be 
extremely useful to study the temperature dependence of these chiral symmetry break- 
ing operators on the lattice and study its phenomenological consequences using these 
constraints. 

3. QCD SUM RULE RESULT FOR ImIT y 

We will now discuss specific changes of Imll y in the QCD sum rule approach. 

3.1. Borel sum rule for II at finite density 

After the Borel transformation with respect to u, the dispersion relation in Eq.(9) with 
the OPE in Eq.(|TTj) looks as follows. 



J [ImllV) " ImnV WH]e-^ M W = £ ( n/2 ), M n & ° ( 13 ) 

Motivated by its form in the vacuum, the spectral density will be modeled with a pole 
and a continuum. 

FT 

Imn°(w) = — — - + C0(U 2 - S ) + Pscattering (14) 

[U Z — ITLyY + (myl Y 

The Borel transformation will reduce the uncertainty associated with approximating the 
continuum with a sharp step function. Nevertheless, it gives a maximum Borel Mass 
M^ ax above which the continuum contribution becomes too large. p sca ttering is the possible 
subtraction constants proportional to S(u 2 ), which can be calculated using the nucleon- 
hole contribution! 28, p9| . 



The next step is to calculate the changes in the condensate. We will use the linear 
density approximation. In this limit, the nuclear matter expectation values of any local 
operator O can be approximated as , 

(0) n . m . = (O) + n n {O) N , (15) 

where (-)o is the vacuum expectation value, (-)jv is the nucleon expectation value and n n 
is the baryon density. Up to dimension 6, all the relevant expectation values appearing in 
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the II sum rules are known except for the dim. 6 four quark operator, for which we will 
use the mean field approximation p8| , p9| . The truncation at dim 6 operator determines the 



minimum Borel mass M^ ax , which is determined by assuming that the power corrections 
are less than 30% of the perturbative contribution. 

The final step is to allow the phenomenological parameters to change also in density 
and determine its change by making a best fit of the left hand side to the right hand side 
of Eq.(13) between the Borel interval M^ in and M^ ax . 

This was first applied to the vector mesons by Hatsuda and Lee [^] in the limit r — > 0. 
The vector meson peak was shown to shift as follows, 

^4 = 1- (0-16 ±0.06)^, 
m v (0) n 

= 1- (0.03 ±0.015)^, (16) 

where the first line denotes the result for the p and u meson and the second line that for 
the 4> meson. n denotes the nuclear saturation density. The result for the p, uj meson 
were due mainly to the changes in the four quark condensate ((qq) 2 ) and about 30% from 
{q^^D v q). The former in the mean field approximation is a chiral symmetry breaking 
operator and the latter is a chirally symmetric operator. As for the <fi meson, the effect 
comes mainly from the chiral symmetry breaking operator (m s ss). Therefore, the shift 
in masses come mainly from chiral symmetry restoration. 



Recently, the work has been extended to include possible change in the width |30|| . It 
was found that the downward shift of the mass could be compensated by an increase in 
the width. Hence, there is a band of region in the mass vs. width plane where the sum 
rule in Eq.(13) is reasonably satisfied. However, there exist a best fit, which starting from 
the vacuum value changes as follows for the p, uj meson, 

Amy 50 — MeV , Ar ~ +200 — MeV. (17) 

n n 

Similar analysis are in progress for the <fi meson. 

Summarizing the results on n°, QCD sum rule suggests that due to chiral symmetry 
restoration, there will be a decrease in mass and increase in the width. 

3.2. Borel sum rule for n 1 at finite density 

As shown in Eq.(8), the dilepton spectrum will have contribution from both the longi- 
tudinal and transverse polarization function, which will now be a function of both u and 
q. We will allow the parameters in Eq.(14) to vary non-trivially by a term proportional 
to q 2 and implicitly the nuclear density n n . Such as, 

F^F + f-q 2 , my + a-q 2 , S -> S + s ■ q 2 . (18) 
In the limit r — > 0, this will give the following form for IulII 1 , 

ImnV) = / • 5(u 2 -ml) -a- F5'(u 2 - m 2 v ) - s ■ c 5(u 2 - S ) + M 2 n n b scatt 5\u 2 ). (19) 

Here, we have included the subtraction constant associated with the S'(u 2 ) singularity, 
which can be calculated from the nucleon-hole contribution. There is another unknown 
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S(u 2 ) singularity. To eliminate this, we use a once subtracted dispersion relation for IT 1 . 

(_l)n/2 



f pMI^u) - Imn 1 (u) pert («)]w 2 e- u2 / M2 (iu 2 = £ 



n=2,4 



(n/2 - 1)!M"+ 2 



hi 



(20) 



From a best fit analysis, one can extract the momentum dependent parameters |27] . The 
results are shown in Table 4. 



Table 4 

Results for the parameters at nuclear matter density. The values are from best fit of the 
Borel sum rule in eq.(20). 





a 


/ 


s 


Transverse p 


-0.108 


0.190 


-0.028 


Transverse u 


-0.081 


0.171 


-0.010 


Longitudinal p, u 


0.023 


0.066 


0.029 


Transverse <fi 


0.004 


0.010 


0.009 


Longitudinal <p 


0.009 


-0.001 


0.009 



As discussed before, a non-vanishing a will shift the average peak position by AM = 
my + aq 2 — my, even if there is no change in the scalar mass my. With the values of 
a obtained, we find that even at q = .5GeV/c, above which point our formalism breaks 
down, the shifts are less than 2% (0.05 %) at nuclear matter density for the p, uj {(f)). This 



is smaller than the expected scalar mass shift of the p, to (~ 20%) and <p (~3%) |28[ and 
justifies (to a first approximation) neglecting the three momentum dependence and the 
polarization effect when implementing the universal scaling laws (Brown-Rho scaling) of 
the vector mesons in understanding the dilepton spectrum in A-A and p-A reaction |31[|. 

The contribution of the longitudinal and transverse polarization to the dilepton spec- 
trum, depends on the angle between the sum and difference of the three momentum of 
the out going dileptons[0. However, after averaging, the contribution of the transverse 
polarization becomes twice that of the longitudinal polarization, as can be seen from 
Eq.(8). Hence, to a good approximation, one can implement the finite q effect into model 
calculations by including only the transverse dispersion relation. Making a linear fit at 
the nuclear matter density, one can parameterize the vector meson mass in medium as 
follows, 

"' / ' ( "" ) 1 (0.023 : 0.007) ( q ^ "' 



m p (0) VO.5/ n 

mM - I- (0.016 ± 0.005) 



77^(0) VO.5/ no 

m ^ nn) - 1 + (0.0005 ± 0.0002) (-^] 2 —, (21) 



m^(0) VO.5/ n 

where q is in the GeV unit. In a full model calculation, the effect in Eq.(21) should be 
added to that of Eq.(16). It should be noted that the main operators responsible for the 
result in Table 4 come mainly from chirally symmetric twist-2 operators. So the dispersion 
effect is a chirally symmetric effect and has little relation to chiral symmetry restoration. 
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4. SUMMARY 

1. Chiral symmetry is restored at high temperature and/or density, while Ua(1) break- 
ing effect persist above the transition. It would be useful to perform lattice calcula- 
tions of C/a(1) breaking effect above the chiral phase transition using 2 quark-bilinear 
order parameter for quark flavors larger than 2. 

2. Using the QCD constraints one can relate observed changes in spectral density to 
changes of local operators at finite temperature or density. A systematic lattice 
gauge theory calculations on the temperature dependence of the chiral symmetry 
breaking condensate would be a useful future project. 

3. QCD sum rules in matter suggests that the vector meson width will increase and its 
peak position shift downwards in energy at finite density and temperature before it 
becomes part of the continuum above the critical boundary. These changes are due 
mainly to chiral symmetry restoration. 

4. QCD sum rule in matter suggests that the dispersion effect ( finite q effect) is small, 
but should be included in a full study of the changes in the vector meson peak in 
medium. 
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